Mon. Not. R. Astron. Soc. 000, [T}{4](2006) Printed 5 February 2008 (MN WT^X. style file v2.2) 



Fuelling Active Galactic Nuclei 



A. R. King 1 and J.E. Pringle 1,2 

1 Theoretical Astrophysics Group, University of Leicester, Leicester LEI 7RH 

2 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CBS OH A 



5 February 2008 



ABSTRACT 

We suggest that most nearby active galactic nuclei are fed by a series of small- 
scale, randomly-oriented accretion events. Outside a certain radius these events pro- 
mote rapid star formation, while within it they fuel the supermassive black hole. We 
show that the events have a characteristic time evolution. This picture agrees with 
several observational facts. The expected luminosity function is broadly in agreement 
with that observed for moderate-mass black holes. The spin of the black hole is low, 
and aligns with the inner disc in each individual feeding event. This implies radio jets 
aligned with the axis of the obscuring torus, and uncorrelated with the large-scale 
structure of the host galaxy. The ring of young stars observed about the Galactic 
Centre are close to where our picture predicts that star formation should occur. 
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1 INTRODUCTION 

' It is now generally accepted that the nuclei of most galaxies 
contain supermassive black holes. Assembling the observed 
hole masses by accretion at plausible efficiency accounts for 

' the emission of active galactic nuclei (AGN) over cosmic 
time (Soltan, 1982; Yu & Tremaine, 2002). This implies that 
essentially every galaxy is intermittently active, in phases 
when its black hole grows. Only a small fraction of galaxies 
are observed to be active (e.g. Heckman et al., 2004), so we 
know that these phases are short. 

Yet we still do not understand how gas gets down to 
the black hole. There are a number of candidates, includ- 
ing galaxy mergers (major and minor), bars, bars within 
bars, turbulence in the ISM, stellar mass loss, and viscous 
accretion discs. However, the nature of the problem is of- 
ten underestimated. For example, in a recent review article, 
Wada (2004) discusses 'Fuelling Gas to the Central Region 
of Galaxies' in terms of bringing gas inwards to around 100 
pc from the central black hole. But, as emphasised by Shlos- 
man, Begelman & Frank (1990), the gas must get to remark- 
ably small radii before the inner viscous accretion disc can 
take over and bring the gas to the hole within the observed 
activity timescale. The inflow timescale through a disc of size 
1 pc is already ~ 10 9 yr, even with maximal assumptions 
about viscosity. For shorter AGN phases, or less efficient vis- 
cosity, the gas must be fed to the central disc at still smaller 
radii. 

Observations show a wide range of activity, from highly 
luminous quasars at redshifts up to z ~ 6, through the 
barely detectable nuclei in LINERS, to the weak activity 



of our own Galactic Centre. It seems likely that all of the 
fuelling processes discussed so far may play a role in at least 
some galaxies. In this paper we focus attention on the fu- 
elling of the nucleus in a normal Seyfert galaxy, with a cen- 
tral black hole of mass M ~ 10 7 — 10 s M , (Ferrarese et 
al., 2001; Heckman et al., 2004; Denney et al., 2006) and 
luminosity L ~ 0.1L_e (with Le the Eddington luminos- 
ity) thus accreting at a rate M ~ 0.2Mq yr -1 . In a typical 
Seyfert event lasting, say, 10 6 yr, this would imply an ac- 
creted mass AM ~ 2 x 10 5 M (e.g. Emsellem, 2002, 2004; 
Martini, 2004) . This is ridiculously small compared with the 
gas available in a typical spiral galaxy, and we immediately 
see that we are not concerned with major gas flows involv- 
ing spiral arms or galactic bars. Feeding the black hole at 
such modest rates needs only some small-scale event to drop 
a little material of low angular momentum into the central 
regions of the galaxy. 

We have two pieces of evidence suggesting what such an 
event might be. First, Kinney et al. (2000; see also Nagar & 
Wilson, 1999) found that in a sample of nearby [z < 0.031) 
Seyfert galaxies, the direction of the jet from the central 
black hole, and therefore presumably the orientation of the 
central regions of the accretion disc, were unrelated to the 
orientation of the disc of the host galaxy. Second, Schmitt 
et al (2003) surveyed extended [O III] emission in a sample 
of nearby Seyfert galaxies, and found that although the [O 
III] emission is well aligned with the radio, there is no cor- 
relation between its orientation and the major axis of the 
host galaxy. Assuming that the orientation of the [O III] 
emission is governed by the geometry of the inner torus, of 
typical radius 0.1 - 1.0 pc (e.g. Antonucci, 1993), this means 
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that the central disc flow has angular momentum unrelated 
to that of most of the gas in the host galaxy. 

This tells us two things (see also the discussion in Kin- 
ney et al., 2000). First, in line with our inference above, if the 
galactic disc supplied the torus material via bar-driven disc 
evolution, or grand-design nuclear spirals, or similar mech- 
anisms, some means of randomising its rotation axis in the 
central regions is needed. For without some such mechanism 
we expect such processes to leave evidence of the angular 
momentum of the galactic disc, which is not seen. There is 
currently no suggestion as to what this mechanism might 
be, although one can appeal to the fact that the galactic 
disc is much thicker than the radial scale of a few parsecs 
that we are interested in. In the absence of a randomising 
mechanism the gas must come from outside the galaxy as 
part of a succession of (very) minor mergers. Where in a 
galaxy a very small merging satellite deposits its gas is not 
straightforward to compute (Velasquez & White, 1999; Tay- 
lor & Babul, 2001). Kendall, Magorrian & Pringle (2003) 
find that if such small merging satellites are to deposit gas 
near the nucleus then their initial orbits must be fairly ac- 
curate shots. Kendall et al. (2003) also conclude that such 
small merging satellites which do manage to reach the nu- 
cleus arrive there on more or less randomly oriented orbits. 

Second, the inner disc (specifying the radio jet axis) 
remains aligned with the torus. In an earlier paper (King 
& Pringle, 2006, Section 3) we showed that repeated small- 
scale fuelling has precisely this effect. It causes counteralign- 
ment of the hole (King et al., 2005) in about one-half of all 
accretion events, and thus spindown, which is more efficient 
than spinup. Once the hole's mass has doubled, its angular 
momentum is always smaller in magnitude than that of any 
accretion event. It then always aligns its spin with the inner 
disc and thus the torus. 

Thus we have argued that the fuelling of low luminosity 
AGN proceeds via a series of randomly oriented, small-scale 
accretion events, acting directly in the region of the central 
black hole. The material deposited in such an event is likely 
to settle quickly into a ring or disc of material within a few 
local dynamical (or orbital) timescales idyn, where 
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It is well known that such a disc is likely to be self- 
gravitating outside some radius R sg ~ 0.01 — 0.1 pc (Shlos- 
man et al., 1990; Collin-Souffrin & Dumont, 1990, Hure et 
al., 1994). Further, in the outer regions of these discs the 
cooling timescales are sufficiently short that self-gravity is 
likely to cause star formation, rather than enhanced angular 
momentum transport (Shlosman & Begelman, 1989; Collin 
& Zahn, 1999). For the discs we consider in this paper, which 
are relatively thin (because of efficient cooling) and of low 
mass (Mdisc <C Afhoic) self-gravity appears first in modes 
with azimuthal wavenumber m ~ R/H, producing transient 
spiral waves which transport angular momentum (Anthony 
& Carlberg, 1988; Lodato & Rice, 2004, 2005). It is reason- 
able to suppose that where the disc is locally gravitationally 
unstable in this way, most of the gas forms stars. 

We therefore propose the hypothesis that all, or at least 
most of, the gas initially at radii R > R sg is either turned 
into stars, or expelled by those stars which do form, on a 



rapid (almost dynamical) timescale (see also Shlosman & 
Begelman, 1989). This corresponds to a nuclear starburst 
of the kind often associated with AGN (e.g. Scoville 2002). 
A starburst like this lasts around 3 x 10 6 yr in terms of 
its ionising flux (O and early B stars). We suppose further 
that all the gas which is initially at radii R < R Bg forms a 
standard accretion disc, which slowly drains on to the black 
hole and powers the AGN. We stress that under the simple 
form of this hypothesis the nuclear starburst and the AGN 
are two different manifestations of the same accretion event, 
but do not feed one another. We note that such a picture 
provides at least a qualitative explanation for the existence 
of the ring(s) of young stars seen around the black hole in 
the centre of the Milky Way (Genzel et al., 2003). 

In Section 2, we investigate the properties of the accre- 
tion disc at radii R < R Bg . In Section 3, we consider the 
time-dependence of such an event and in Section 4 present 
a simplified luminosity function under the assumption that 
all events are identical, but randomly timed. 



2 PROPERTIES OF THE DISC 

We use the steady disc properties derived by Collin-Souffrin 
& Dumont (1990) in the context of AGN, which are es- 
sentially the same as those derived by Shakura & Sunyaev 
(1973) for steady discs in the context of X-ray binaries. The 
disc surface density is given as 
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Here a is the standard Shakura & Sunyaev (1973) viscosity 
parameter, e is the accretion efficiency, so that the luminos- 
ity L and accetion rate M are related by 



L = eMc 2 , 
also 

L E = 1.4 x 10 46 M 8 ergs _1 , 



(3) 
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is the Eddington luminosity, Ms is the black hole mass, M, 
in units of 10 M B , R is the radius and 



R s = 2.96 x 10 M 8 cm 



(5) 



is the Schwarzschild radius of the central black hole. 

Then the mass of the disc M(< R) interior to radius R 
is given by 
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The disc semi-thickness H is given by 

r = L94 x 10 (m) (al) 
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The condition for the disc to become self-gravitating is 
approximately (e.g. Pringle, 1981) 



(7) 
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We note that this implies 
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almost independently of the black hole mass. This arises 
because (|10|) is an integrated form of the standard steady- 
state disc equation M — 3nvY, combined with (|8j) , which we 
can express as 
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where c s is a mean sound speed and we note from (0 that 
H/R ~ constant. Encouragingly, we see that 7? sg is only 
slightly smaller than the inner edge R ~ 0.03 pc of the ring 
of young stars seen around the black hole in the centre of 
the Milky Way (Genzel et al., 2003). This is to be expected, 
as the disc within _R sg must pass its angular momentum to 
the self-gravitating region further out which in our picture 
produces these stars. 

The mass inside the radius i? sg is given by 
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which is of course just (H/R)M with H/R evaluated at i? sg 
(cf. equation [JJ. 

The accretion rate is given by M — L/ec 2 , which implies 
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Then the evolution timescale of the disc is given by r sg 
M sg /M, which gives 
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Again we note that r sg ~ (H/R)(M/M). 



3 PROPERTIES OF AN ACCRETION EVENT 

In the previous Section we established the general properties 
of the non-self-gravitating disc, assuming that it was ac- 
creting steadily, and giving rise to a luminosity L. Of course 
this is not precisely what we require for the present prob- 
lem. What we are assuming is that the disc, of initial mass 
M sg , is deposited by the accretion event in some unknown 
configuration, and then evolves due to viscosity. As for most 
accretion discs it is fair to assume that initially most of the 
mass, and most of the angular momentum of the disc, is 
predominantly at large radius, and therefore at around i? sg . 
But whatever the initial configuration, within a timescale of 
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Figure 1. Form of the AGN luminosity function predicted by 
the fuelling process, and subsequent disc evolution, discussed in 
this paper. The fraction F, of those sources with luminosities less 
than L/Le is shown as a function of L/L^. This is similar in 
form to those presented in Figure 3 of Heckman et al., 2004 



at most r sg given by ea. (|14p the disc evolves to resemble a 
steady disc at radii R < R sg (e.g. Pringle, 1981). Thus using 
the formulae derived in Section 2 gives a good approxima- 
tion to the actual disc properties after a brief initial period. 
We can now use this to estimate the properties of the disc, 
and therefore of the AGN event, as it evolves. 

From eqn.© we see that E tx M 3/5 R~' i/5 . For a steady 
disc we know (Pringle 1981) that M tx z^E, where v is the 
viscosity. From these two relationships we can deduce that 
for these discs v tx Y?^R. For such discs, similarity solutions 
(Pringle 1991) imply that at late times the accretion rate, 
and hence luminosity, varies with time as L oc t -19//16 . Thus 
we expect luminosity evolution roughly to follow 
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4 IMPLICATIONS FOR THE LUMINOSITY 
FUNCTION 

Heckman et al. (2004) find that most of the current accretion 
on to black holes is on to those with masses in the range 
10 r — 10 8 Mq. In addition they find that of low mass black 
holes, with masses M < 3 x 10 7 Mq, only 0.2 per cent 
are growing at a rate which accounts for 50 per cent of the 
fuelling. They point out that this implies that strong fuelling 
therefore only lasts for a time tf ue \ ~ 0.002 x tn, where 
tu ~ 1.4 x 10 10 years is the age of the universe. Thus the 
total time during which strong fuelling is taking place is 
around tf tte ] ~ 3 x 10 7 years, or, if there are fuelling events 
per black hole per Hubble time, it implies that each event 
lasts around t av ~ 3 x 10 7 /N yr (Heckman et al., 2004). 
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Accordingly we make some simplifying assumptions to 
work out the luminosity function we might expect from our 
model. We take a = 0.03 and e = 0.1, and consider black 
holes with mass 1O 7 M as representative of the range mak- 
ing the major contribution to the distribution of accreting 
black holes. We assume that these black holes undergo ran- 
dom, but identical, accretion events such that the initial 
luminosity of each event is Le- In this case, we see from 
the above that the initial disc mass is M sg ~ 5.95 x 10 4 
Mq, and the initial evolution timescale is r sg ~ 2.41 x 10 5 
yr. Equating this timescale to the length of each event, 
i C v we see that the number of events has to be about 
TV = 0. 002t£r/T S g ~ 116, and therefore the average time 
between events is t rcp ~ tn/N = r sg /0.002 ~ 1.2 x 10 s yr. 

Writing / = L/Le, our assumption implies that the 
initial value of / is 1. The average final value is / cn d = 
/in/(l + A"^ 16 ), where A cnd = t rop /r sg = 500, and so the 
maximum possible range of / is (1 + A ell d) 19//16 . Of course 
observations cannot probe this full range, and from Figure 
3 of Heckman et al. (2004) we see that for 10 7 M black 
holes the observed range of / is around 40 (corresponding 
to A > 18 and N < te/18r sg = 540). 

Inverting equation \T5\ we find that 

— = r l6/l9 - 1, (i6) 

for < t < r sg . Then assuming that fuelling events for dif- 
ferent black holes are independent, we find that the fraction 
F(> /) of black holes with luminosities > / is given by 

r-16/19 _ -, 

F (> /) = ;_ 16/19 • (i7) 

Tend 1 

We plot this in Figure 1. As can be seen, this curve has 
a similar form to the distributions of low-mass black holes 
(3 x 1O 6 M - 3 x 1O 7 M ) in Figure 3 (Left) of Heckman 
et al. (2004) in the observed range 1 < / < 40. Given the 
simplicity of our assumptions (identical independent fuelling 
events etc) we regard this as encouraging. 



5 CONCLUSIONS 

We have suggested that the feeding of most nearby ac- 
tive galactic nuclei proceeds via a series of small-scale, 
randomly-oriented accretion events, rather than large-scale 
events bearing the imprint of the host galaxy. Outside a 
certain radius these events cause rapid star formation, while 
within it they feed the supermassive black hole. This picture 
implies a characteristic time decay of each event and implies 
a luminosity function broadly in agreement with that ob- 
served for moderate-mass black holes. The chaotic nature 
of the accretion keeps the black hole spin low, allowing each 
individual feeding event to produce radio jets aligned with 
the axis of the obscuring torus, which itself is uncorrelated 
with the large-scale structure of the host galaxy. Our picture 
predicts that star formation should occur at radii compara- 
ble with those of the ring of young stars observed about the 
Galactic Centre. In an earlier paper (King & Pringle, 2006) 
we showed that small-scale feeding events allow supermas- 
sive black holes to reach large (>1O 9 M ) masses at high 
redshifts ~ 6. We conclude that small-scale chaotic feeding 



offers a promising explanation of the growth of most su- 
permassive black holes, and thus for the activity of galactic 
nuclei. 
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